Paraxial light in a Cole-Cole nonlocal medium: integrable 

regimes and singularities. 

Boris Konopelchenko and Antonio Moro* 

Dipartimento di Fisica dell'Universita di Lecce 
and INFN, Sezione di Lecce, via Arnesano 1-73100 Lecce, Italy 



ABSTRACT 

Nonlocal nonlinear Schrodinger-type equation is derived as a model to describe paraxial light propagation in 
nonlinear media with different 'degrees' of nonlocality. High frequency limit of this equation is studied under 
specific assumptions of Cole-Cole dispersion law and a slow dependence along propagating direction. Phase 
equations are integrable and they correspond to dispersionless limit of Veselov-Novikov hierarchy. Analysis of 
jj^ . compatibility among intensity law (dependence of intensity on the refractive index) and high frequency limit of 
Poynting vector conservation law reveals the existence of singular wavefronts. It is shown that beams features 
depend critically on the orientation properties of quasiconformal mappings of the plane. Another class of wave- 
fronts, whatever is intensity law, is provided by harmonic minimal surfaces. Illustrative example is given by 
helicoid surface. Compatibility with first and third degree nonlocal perturbations and explicit solutions are also 
discussed. 
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> ; 1. INTRODUCTION. 

Propagation of the light through spatially nonlocal media is one of main subjects in modern nonlinear optics. 1-6 
Recently, it has been shown, both theoretically 3 and experimentally, 4, 5 that highly nonlocal Kerr-type media 

\& supports propagation of stable laser beams, usually referred to as "solitons" . Laser beams propagating in Kerr- 
type media are modelled by the celebrated (2 + 1)D nonlinear Schrodinger (NLS) equation. Unlike (1 + 1)D 
NLS equation, it is not integrable by inverse scattering method. 7 In the case of high nonlocal media, (2+l)D 
NLS equation can be approximated by quantum harmonic oscillator (QHO) equation. Using results concerning 
QHO in the optics context, it has been possible to show the existence of solitons. 3 Nevertheless, a complete 

r— h : understanding of the role played by spatial nonlocal effects in the propagation of radiation in nonlinear media 
is still an open problem. Hence, it is of interest to investigate possible integrable regimes for three-dimensional 
nonlinear models. 

^ ■ In the present paper, we propose the nonlocal nonlinear Schrodinger-type (NNLS) equation as the model to 

describe a beam propagating in a medium with a certain spatial nonlocality degree and obeying some intensity 
law. By "intensity law" we mean certain functional dependence among dielectric function and intensity of the 
beam. We study the high frequency limit under specific assumptions, namely, the Cole-Cole dispersion law for 
the medium and suitably slow variation of the fields along propagating direction. Separating behaviors on the 
transverse and propagating directions, we derive two couples of equations for the phase and the intensity. Due 
to the intensity law, phase and intensity equations are strongly coupled and their compatibility appears to be a 
highly nontrivial problem. 

As far as phase equations concerned, the compatibility has been studied in the papers [8,9]. It has been 
shown that compatibility condition corresponds to the request that refractive index satisfies an infinite set of 
integrable equations referred to as dispersionless Veselov-Novikov (dVN) hierarchy. In particular, degrees of 
spatial nonlocality are in one-to-one correspondence with the orders of equations of dVN hierarchy 
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So, we study compatibility among phase and intensity equations for certain assigned intensity law. We focus 
first on transverse equations. Their compatibility leads us to the study of a second order nonlinear partial 
differential equation. Since for most of physical models it satisfies the ellipticity condition we restrict ourselves 
to the class of so-called elliptic intensity laws. 

The use of known results in the theory of elliptic equations and the special properties of Beltrami equation 
allows us to conclude that, for a generic elliptic intensity law, there exists singular wavefronts which can be 
used to describe self-guided beams. Moreover, we demonstrate that harmonic minimal surfaces provide us with 
a class of wavefronts for any elliptic intensity law. As explicit examples, helicoidal wavefronts are discussed. 
We stress that the helicoid is a well known example of singular wavefront which possesses several interesting 
properties. 10-16 We show that helicoidal wavefronts are preserved for a specific class of first degree nonlocal 
perturbations up to next-to-leading order contributions. We present a simple example where the first degree 
nonlocality is responsible of a slight stretching or compression of the hclicoid's pitch. 

General problem of compatibility among equations for transverse and propagating direction is difficult as 
much as crucial. In our model the consideration of propagating direction contributions is equivalent to inclusion 
of nonlocal terms. Besides, in the local case, equations for propagating direction are automatically compatible 
for any intensity law and it is sufficient to consider the transverse equations. We demonstrate that there exists a 
class of solutions for the third degree nonlocality which are in agreement with a suitable intensity law. Exploiting 
hydrodynamic-type reductions, found in the paper [17], we present some explicit formulae. 

The paper is organized as follows: in section 2 derivation of nonlocal nonlinear Schrodingcr equation is 
presented. In the section 3 phase equations and high frequency limit of Poynting vector conservation law in a 
Cole-Cole medium are displayed. A study of compatibility among intensity law and transverse equations and 
properties of solutions is presented in the section 4. Compatibility of intensity law and nonlocal perturbations 
is illustrated in section 5 by an explicit example. Some concluding remarks close the paper. 

2. NONLOCAL NLS-TYPE EQUATION. 

Let us consider the nonlinear Maxwell equations 

<9B 

V A E+ — = V • D = 

dt (1) 

V A B — =0 V-B = 

at 

where the displacement vector D = D(E) depends nonlinearly on the electric field. We look for time-oscillating 
solutions such that 

E (x, t) = E (x) e~ iajt , D (x, t) = D (x) e~ iajt . 
Under these assumptions one gets the following second order equation 18 

V 2 E + w 2 D- V(V-E) = 0. (2) 

Displacement vector D is chosen in a way to collect both nonlinear and nonlocal effects, i.e. 

D = eE + Fnloc- (3) 

First term of the r.h.s. contains the nonlinearity. Dielectric function e is assumed to contain a weakly nonlinear 
perturbation of the form 

e = e a + a 2 e, a « 1, (4) 

where e depends on the intensity of the radiation. 

The second term in the r.h.s. of (3) provides us with the spatial nonlocal response of the medium. We assume 
it to be of the form 

n— 1 /,m— 1 / 
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One gets the expression (5) for the nonlocal term considering an integral relation among the displacement and the 
electric field where the kernel is a sum of 5-functions and spatial derivatives. We define the degree of nonlocality 
as the maximum derivative order which appears in expression (5). 

Paraxial approximation is performed according to the rule (see e.g. Ref. 19) 

E -» aEe ikz D -» aBe ikz (6) 

(setting vacuum light speed c = 1, k = ^/iou), and the amplitudes are assumed to depend slowly on the 
coordinates in such a way that the gradient acts as follows 



dx' dy 1 dz J \ ' dz 

Using the expression (3) and the rules (6) and (7) in the second order equation (2), one shows that the first 
nontrivial contribution is of the order O (a 3 ) and it provides us with the nonlocal nonlinear Schrodingcr-type 
(NNLS) equation 

<9E 

2i^u— + ViE + w 2 eE + T NLoc = (8) 

where 

Fnloc = Tnloc + O (a 4 ) . (9) 

Note that in absence of nonlocal contribution Tnloc and for a Kerr type nonlinear response of the medium, 
equation (8) becomes the standard nonlinear Schrodinger equation in 2 + 1— dimensions. 

3. GEOMETRICAL OPTICS LIMIT OF NNLS EQUATION IN COLE-COLE MEDIA. 

In this section we study the geometrical optics limit of NNLS equation for a medium which satisfies so-called 
Cole-Cole dispersion law 20 

6 1 

1 + {itqu) zv 2 

where To is the relaxation time and v is a suitable phenomcnological exponent whose range of values is of crucial 
importance in the following derivations. Relation (10) describes frequency dependence of dielectric function for 
several polar media as well as liquid crystals and various biological tissues (see e.g. Refs. [20,21]). 

We study a particular high frequency regime of NNLS equation where the phase of electric field varies slowly 
along z-direction in such a way that 

We also assume nonlocal effects to be weak in high frequency limit in such a way that they contribute at 
most with uo~ 2v — order terms 

Tnloc = u- 2v f NLoc + O (12) 
In order to calculate geometrical optics limit we adopt, as usual, the following representation of electric field 

E = f(x,y,cj 2v ) e^Ow 2 "*). (13) 

Expansion of all functions with respect to small parameter lo~ 21 ' 

f(x, y, uj 2u z) = i(x, y, z) + w'^fi (x, y,z) + ... 

S(x, y, u? v z) = S(x, y, z) + w~ 2v S 1 (x, y,z) + ... (14) 
i(x, y, u 2v z) = i(x, y, z) + uj~ 2ij e x (x, y,z) + ... 
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together with dispersion law (10) and the prescriptions (11) leads us to the high frequency limit of NNLS equation 



,2-2;/. 



(V ± 5) 2 = 4u, (15a) 
S z = tp(S x ,S v ,x,y,z) (15b) 

where Au = e and ip = (ii + TnlocJ / (2y / £o). For convenience we denote the partial derivative of a function 

/ with respect to variable i] as f n . We stress that condition on the exponent < v < 1/2 is crucial to 
separate equation (15b) from polarization contributions. Equation (15a) is the standard eikonal equation in 
two-dimensions while the function ip in equation (15b) is an N-degree polynomial in S x and S y and it describes 
an N-degree nonlocal response. Moreover, in agreement with inversion phase symmetry (S — > —5) of the eikonal 
equation, polynomial degree of function ip must be odd. We note also that the system of equations (15) has been 
derived first directly from the Maxwell equations. 9 In the present paper, we will consider local response and 
first and third degree nonlocal responses. 

In local case we have 

<p = tp(z), u = u(x,y), (16) 
where ip and u are certain function of their arguments. For a first degree nonlocality we have 

p = a 1 S x + a 2 S y , (17) 

and u must satisfy the following linear equation 

u z = (a lU ) x + (a 2 u) y (18) 

where at\ and a 2 are harmonic conjugate functions (they satisfy Cauchy-Riemann conditions a.\ x = a 2y and 
ai y = —U2x)- Finally, for a third degree nonlocal response we have 

P=\sl-\s x S 2 y + V 1 S x + V 2 S y (19) 

and u satisfies the dispersionless Veselov-Novikov equation 

u z = (V lU ) x + (V 2 u) , 

V lx - V 2y = -3u x , (20) 

Vly + V 2x = 3Uy. 

In order to complete geometrical optics description we construct high frequency limit of Poynting vector 
conservation law 

V-P = 0. (21) 



Poynting vector can be written down as follows 



VS 1 * 

P = I—— (22) 



where 

S* = ^ z + S{x,y,uj 2v z) (23) 

is the total phase of electric field and I is the intensity. In the following, for sake of simplicity, we set dielectric 
constant eq = \. Using the previous prescriptions to perform high frequency limit in paraxial approximation and 
expanding intensity 

I = I +LO- 2v h, (24) 
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one concludes that equation (21) leads to the following set of equations for I and I\ 

V ± / o -V ± S' + /oViS' = (25a) 

V±Ji • V±S + I 1 V 2 ± S+- S ± = 0. (25b) 

oz 

In this regime, properties of the light are completely described by the systems of equations (15) and (25). If 
no additional conditions are considered, once compatibility of equations (15) is guaranteed, linear system (25) 
provides us with leading orders intensity contributions I and I\. 

Integrability of the phase equations (15) has been discussed in details in Refs. [8,9]. In particular, it has 
been shown that for a function <p with a polynomial dependence on S x and S y , compatibility condition of 
equations (15) selects the class of refractive indices which satisfy the dispersionless limit of the Veselov-Novikov 
(dVN) hierarchy. In particular, dVN hierarchy allows us to classify different degrees of nonlocal response by a 
one-to-one correspondence among degrees of the polynomials ip and the equations of the hierarchy. 

We mention that dVN hierarchy is amenable by the quasiclassical <9-dressing method 9 and a reduction method 
based on symmetry constraints appears to be an effective method to construct exact solutions. 17 



4. PHASE AND INTENSITY TRANSVERSE EQUATIONS. 
4.1. Elliptic intensity law. 

Several models in nonlinear optics analyze a self-action of the medium assuming that the influence of an external 
electromagnetic field is described by a functional dependence among intensity and refractive index I = /(it). This 
is the form of intensity law in the geometrical optics limit. Kerr type media, where / oc u, and saturable nonlinear 
media such that u oc log(l + I/It), where I t is the so-called threshold intensity, are two famous examples. Of 
course, for self-action nonlinear media models phase equations (15) and intensity equations (25) appear to be 
strongly coupled and the study of their compatibility is a highly non-trivial problem. In the present section, we 
analyze compatibility and properties of equations (15a) and (25a), here referred to as transverse equations. We 
note that this analysis is quite trivial in the case of the media exhibiting local response. Indeed, according to 
relations (16) u does not depend on the z-variable. Then, equation (25b) coincides with equation (25a) and they 
are automatically compatible with (15b). 

Let us consider a self-action nonlinear medium obeying intensity law Iq = Iq(u). Observing that V±Iq = 
IqV±u, where "prime" means the derivative with respect its argument u, and using equation (15a) in (25a), one 
gets the following second order partial differential equation 

AS XX + BS yy + 2CS xy = 0, (26) 

where by definition 1 = log Iq and 

A = l 'S 2 x + 2, B=l 'S 2 + 2, C=1 'S X S V . (27) 

Recall that properties of a second order partial differential equation depend critically on the discriminant 

A = AB-C 2 . (28) 

If A > 0, equation is of the elliptic type when A > , parabolic when A = and hyperbolic when A < 0. In 
this paper we focus on the elliptic case motivated by the fact that a wide class of such models for optical spatial 
solitons is provided by a Kerr dependence or logarithmic saturable nonlincaritics. Indeed, one can simply verify 
that for a Kerr-type intensity law of the form I = vr where j3 > and a logarithmic one, u = log(l + /o//o,t), one 
has A = 4(22o 'u + 1) > uniformly. Ellipticity of equation (26) seems to be associated to important physical 
properties of a wide class of nonlinear media. Moreover, elliptic equations such as (26) have been well studied 
in mathematics and possess several remarkable analytical and geometrical properties. 22-27 
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Now, we will study equation (26). Introducing complex variable A = x + iy and the complex gradient 
w = S x — iS y , one can rewrite equation (26) cquivalently as follows 



aw\ + bw^ + awi + bw\ = 



(29) 



where 



a 



= \{A-B + 2iC), b=\ 



(A + B — 2iC) . 



(30) 



A class of solutions of equation (29) is provided by the solutions of, the so-called, nonlinear Beltrami equation 



It is straightforward to verify that ellipticity condition A > implies < 1. Beltrami equation describes 
so-called quasiconformal mapping of the plane of complex dilatation /z. 26,27 This fact highlights an intriguing 
connection among geometrical optics and quasiconformal mapping theory. 

For instance, let us consider a laser beam whose input profile is Iq ^ inside a simply connected domain G 
and Io = on the smooth boundary r of G and outside it. Because of functional relation Iq = Io(u), refractive 
index is constant on T too. More specifically, we have 



Boundary T is mapped on the 2y / uo-radius circle. In virtue of equation (32), for A G T one can write w = 
2y / uo exp(i0). Then, if w is assigned in such a way that, for instance, w(0) = 0, where e G is assumed, 
variation of argument of complex number w is 



Under mentioned conditions, it can be demonstrated that w is a homeomorphic mapping of domain G onto 
2y/uQ— radius disk T '. 23 As a consequence, mapping w preserve the topology of domain G. It could be 
interesting to note that if w = w(X, z) "evolves" along z according to equation (15b) and this evolution is 
compatible with equation (26), inverse mapping (it exists since w is homeomorphic) A = \(w), with w £ T ', 
provides us with the beam profile at any z. 

Mapping w — S x — iSy — {Sx-, — Sy) can be also regarded as a two-dimensional vector field on the A — plane 
associated with transverse components of the wavefront normal unit vector. In the case ip = 0, it is has the 
simple form 



From this point of view, mapping w gives us information about light-rays distribution around propagating 
direction. 

As specific example, let us consider a mapping as in figure la. Curves T and T' are oriented leaving domain 
on the right hand side. Under conditions mentioned above, there exists a homeomorphism w of domain G 
onto G' acting in such a way that w(s\,0) = (—1,0), ui(0,s 2 ) = (0,1) and w(— s 3 ,0) — (1,0). Consequently, 
normal unit vector to wave-front on the boundary T points inside the domain G (recall that y-component of w 
reverse y-component of n). This mapping allows the radiation to be 'self-guided' around propagating direction. 
Conversely, if w(si,0) = (1,0), w(0, s 2 ) = (0,-1) and w(— s 3} 0) = (0,-1) radiation spreads transversely far 
from z— axis. Note that in both cases homeomorphism w is sense-reversing. If we consider a sense-preserving 
mapping such as, for instance, w(si,0) = (1,0), w(0,s 2 ) = (0,1) and w(— s 3 ,0) = (—1,0), the beam spreads 
along x— direction and tends to be trapped along y— axis (see figure lb). We emphasize, finally, that all of these 
observations can be generalized to the case of arbitrary n— connected domains. 

Another class of solutions of equation (29) can be obtained solving the equation 




(31) 



w(X)w(\) = 4w, , a e r. 



(32) 



A argw = 2tt. 




(33) 



bw^ + aw^ = 0. 



(34) 
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Figure 1. a) Sense-reversing mappings describe beams which tend to be trapped (solid line) or spread (dashed line), b) 
If the mapping is sense-preserving light beam spreads along a;— direction and tends to be trapped along y— direction. 



Passing to reciprocal coordinates by inversion of the system 



A = X(w, w) 
A = X(w, w) 



(35) 



one converts equation (34) into the form of the linear Beltrami equation 



Au, = v(w, w)X 



(36) 



where v(w,w) — a/b and, due to ellipticity, |^| < 1. The advantage of the use of the reciprocal transformation 
is that the Beltrami equation is just linear and several properties of the analytic functions can be rigorously 
generalized to its solutions (see e.g. Ref. [25]). An interesting result is provided by the generalization of Liouville 
theorem. Indeed, if A = X(w, w) is bounded on whole w— plane and satisfies linear Beltrami equation (36) it can 
be shown (Vekua's theorem) that X(w,w) = constant. 28 Of course, for the constant solution, mapping from 
w— to— A— plane is singular and reciprocal transformation (35) is not defined. Then, any non-trivial solution of 
equation (36) must be singular somewhere on the complex plane and different type of singularities can occur, 
such as poles, essential singularities, singularities of the derivatives etc. A complete analytical and geometrical 
classification of them will be done elsewhere. 29 

For sake of simplicity, let us consider a solution A(w, w) having a simple pole at w = 0. As discussed above, 
one can always consider a homeomorphism from D' = C \ D e to Dr where D e is a disk of arbitrarily small radius 
e on the w— plane and Dr is the i?-radius disk on the A— plane, mapping the boundary of D e on the boundary of 
Dr. Inverse mapping w : Dr , — > D' constructed in such a way can be used to describe a beam "confined" around 
z— axis. Indeed, light rays crossing the boundary of Dr can settled down arbitrarily parallel to z— axis as much 
as disk D e is small. 

Coming back to nonlinear Beltrami equation (31), we expect that for "mild enough" complex dilatations ji, 
Vekua's theorem still holds. In these cases, the only one bounded solution on whole A— plane is w = constant. 
In most of meaningful physical situations, intensity distribution on the A— plane, at certain z, goes to zero for 
A — > oo, or equivalcntly, one can says that intensity vanishes outside a big enough R— radius disk Dr. Thus, 
outside Dr, refractive index assumes a constant value u = uq and the solutions of eikonal equation (15a) is 



where Co, c\ and C3 are constants and the condition Cq + c 2 = Auq holds. For a paraxial beam we have Co = c\ = 0. 
By a consequence w = S x — iS y = c — ic\ = and, of course it satisfies Beltrami equation in C \ Dr. In virtue 
of Vekua's theorem, the only one bounded solution is w = 0. Then, any non-trivial solution must be singular 
somewhere on the plane. In our example, wavefront is approximately plane for A e C \ Dr and possesses 
singularity inside D R . 



S = c x + ciy + c 3 



(37) 



4.2. Minimal surfaces: the helicoid example. 

Standard approach to calculate solutions of Beltrami equation is based on its reduction to a two-dimensional 
integral equation 25, 30 which is amenable by using the successive approximation method. In general, to find 
'explicit' exact solutions is a quite challenging problem and successful chances are critically depending on the 
form of the complex dilatation. Incidentally, we note that solutions possessing cylindrical symmetry such that 

S = S(r), u = u(r), r = ^x 2 +y 2 , (38) 

are not compatible with intensity law. It is straightforward to verify that equations (15a) and (25a) along with 
assumptions (38) imply that intensity depends explicitly on z— axis distance r 

h = T ^- r , (39) 

where c is an arbitrary constant. 

Here, we discuss solution of the equation (26) in connection with the minimal surfaces. They are described 
by the following elliptic equation (see e.g. Ref. 31) 

(1 + Si) S xx + (1 + S 2 X ) S vv - 2S x S y S xy = 0. (40) 

Now, we restrict ourselves to the class of solutions of equation (26) which are also harmonic, that is 

SXX + Syy = 0. (41) 

Using condition (41) in equation (26), one gets the equation 

S X S XX SySyy + 2S x S y S xy — 0, (42) 

for any intensity law. It is straightforward to check that equations (42) and (40) coincide for harmonic solutions. 
In other words, a class of solutions of equation (26), whatever is dependence Iq = Io{u), is just given by the class 
of the harmonic minimal surfaces. An interesting non-trivial example of solution is 



K arctan 



(|), (43, 



where K is an arbitrary constant. It is straightforward to check that function S in (43) satisfies equations (40), 
(41) and (42) simultaneously. Equation of corresponding wavefronts is 

S* = z + K arctan (|) = const, (44) 

where the constant K is the "pitch" of the helicoid. We stress that singular wavefronts of the form (44) are well 
known in the theory of nonlinear singular optics and their topological properties have important phenomcnological 
consequences connected with description of screw wavefront dislocations or so-called optical vortices 12-14 (see 
also Ref. [15,16] and references therein). Note that complex gradient associated with the helicoid 

w = —i^- (45) 
A 

is analytic and it has a simple pole singularity at A = 0. Nevertheless, normal vector to the wavefront, whose 
components coincides (up to a sign) with real and imaginary parts of w, is not defined. Indeed, vector w = 
(S x ,—S y ) has no limit as (x, y) — > 0. Figure 2 shows examples of helicoidal wavefronts (44) usually parametrized 
as follows 

x = vcosi, y — vsint, z = Kt. (46) 
One-started helicoid shown in figure 2a is obtained for < v < +oo; two-started helicoids shown in figure 2b, c 
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b) 



c) 





Figure 2. Helicoidal structures of wavefront around 2— axis: 
(K = 1); c) two-start left screw (K — —1). 



a) one-start right-screw (K = 1); b) two-start right screw 



are obtained for — 00 < v < +00. Refractive index 



K 2 



4(x 2 + y 2 ) 



(47) 



has cylindrical symmetry around z— axis and displays divergence at A = 0. Thus, in most of physical cases 
one has a divergence of intensity Iq. This means that around A = geometrical optics approximation fails and 
wave effects become relevant. In particular, necessary condition for the existence of singular wavefronts is that 
intensity vanishes where phase function is singular. Indeed, in this region interference phenomenon is no more 
negligible and can realize this condition. 

Finally, we note that helicoidal wavefront is 'stable', up to lo~ 2i/ — orders deformations, for a specific class of 
first degree nonlocal perturbations. In this case, the phase (43) is defined up to an additive arbitrary function of 
z— variable i[)(z), such that ip' — <p(z). Thus, using the phase expression S = K arctan (y/x) + ip(z) along with 
relation (47) in equations (17) and (18), we conclude that the helicoid is deformed for nonlocal harmonic data 
a\ and a 2 satisfying simultaneously the following equations 



(x 2 + y 2 ) a lx = xai + ya 2 
—ya x + xa 2 = <p(z) (x 2 + y 2 ) 



(48) 



Trivial solutions a\ = a 2 = = corresponds to the local case discussed above. A simple non-trivial solution 
a\ = —jy, a 2 = "fx and ip = 7, where 7 is an arbitrary constant, provides us with the following wavefront (we 



remind that to return to speed variables one has to substitute ip(z) — > i^(uj 2v z)) 



S* = z + 



K 



1 + u 



-1v 



arctan 
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(f) 



const. 



(49) 



Note that for 7 = equation (49) coincides with the equation (44). As a consequence of nonlocal response the 
helicoid's pitch is compressed if 7 > and stretched if 7 < 0. 



5. INTENSITY LAW AND NONLOCAL PERTURBATIONS. 

Compatibility condition between dVN hierarchy and the intensity law Jo = Jo (u) is an as challenging as intrigu- 
ing issue. Indeed, due to the Beltrami equation, it could establish a direct connection between dispersionless 
integrable systems and an entire class of quasiconformal functions parametrized by z. Unfortunately, intensity 
law appears to be quite restrictive making the problem of compatible nonlocal responses rather non-trivial one. 
This section is devoted to demonstration of the existence of a suitable intensity law such that systems (15) 
and (25) are compatible for a third degree nonlocal response. 
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To do that, we consider hydrodinamic type reductions of dVN equation. They have been found using 
symmetry constraint of the form V^S 1 = u x . 17 It can be shown that dVN equation (20) is reduced to following 
hydrodynamic-type system 



Pi 

Pi 

Pi 
Pi 



1 

2pi-l 2p 2 

An A 12 
A21 A22 



Pi 

P2 



Pi 
Pi 



(50) 
(51) 



where 



An = 3pi (pi - 1) 
A21 = 'ipi (2pi - 1) 



A12 = 3p 2 , 
A22 = 3pi (pi - 1) + 6P2 



and pi := S x , p 2 := 

Now, looking for solutions such that P2 = P2(pi), Pi andp2 are given in implicit form in terms of the following 
algebraic system 



x + G'y + H'z — <& (pi) = 



Pi = 



2c -log (e+v/ITF) 



(52) 
(53) 



£ = p 2 ± dpi + 2 Pl - 1 



where c is an arbitrary constant, G = P2(pi), H = p\ — |pf + |pi(Pi) an d 'prime' means the derivative with 
respect to pi . 

Differentiating eikonal equation (15a) with respect to x and taking into account that u x = V\S, one gets 



~2 



V± ( ) • V ± 5 + Vi5 = 0. 



(54) 



Comparing equation (54) with intensity equation (25a), one obtains the following nice relation among intensity 
and pi— component of the gradient 

Io = Ce~^, (55) 
where C is an arbitrary real constant. Finally, eikonal equation provides us with intensity law 



u(I ) = log 



C 



p 2 -2 log 



(56) 



where last term in r.h.s. is given by algebraic relation (53). 



6. CONCLUSIONS. 

NNLS equation (8) has been proposed as the generalization of NLS equation in order to include different degrees 
of nonlocal responses. Of course, just like (2+l)D NLS equation, NNLS equation is not integrable. Particular 
high frequency regimes discussed in this paper allows us to investigate interplay among a wide class of nonlinear 
responses and nonlocal perturbations free from interference effects. At leading order non-trivial wavefronts ex- 
hibit singular vortex-type behaviors. Moreover, singular phases can support propagation of self-guided beams. 
Harmonic elliptic surfaces provide also with stable beams with singular wavefronts such as the helicoid. More- 
over, there exists special first degree nonlocal perturbations preserving helicoidal wavefronts. In the example 
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considered, first degree nonlocal perturbation produces a stretching or compression of helicoid's pitch. A classi- 
fications of singularities and their properties should provides us with a deeper understanding of this phenomena 
and seems to be a promising subject of study. 

On the other hand it has been shown that intensity law is compatible with third degree nonlocal perturbations. 
In particular, we note that new hydrodynamic-type reductions, obtained by symmetry constraints approach can 
be very useful to calculate new compatible intensity laws. This will stimulates our efforts in that direction. 
Finally, we observe that although the class of solutions in nonlocal perturbation regimes appears to be not 
too large, example discussed above shows they are highly non-trivial. Hence, in this regime, interesting and 
unexpected phenomena could occur. 
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